Charge ordering in oxides : a conundrum solved by resonant diffraction 
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We show that in mixed-valence 3d transition metal oxides undergoing a structural transition, the 
low temperature phase results from an effective ordering of the charge. This arrangement and the 
quantitative evaluation of the atomic charges are determined by using resonant x-ray scattering 
experiments further analyzed with the help of ab initio calculations of the corresponding scattering 
factors. We have found that this reorganization concerns only a small fraction of electron and is 
necessary to reconcile all the experimental data. 
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There is an issue in understanding the charge state 
in the high and low temperature phase of mixed-valence 
transition metal compounds. Few experimental tech- 
niques can access the actual charge state of a given atom 
and most of the accumulated knowledge results from the- 
oretical models pj and from indirect results based on 
transport property measurements Q or from heuristic 
bond- valence calculations on the crystal structure Q. 
Nevertheless, the answer to very simple questions would 
turn out to be of large relevance in the context of metal- 
insulator or insulator-insulator phase transitions issued 
from a charge ordering (CO), chiefly if one could deduce 
a general pattern. For instance, given an average valence 
state as deduced from charge neutrality considerations of 
the chemical composition, would the ground state be an 
effective CO ? and if so how are the corresponding elec- 
trons arranged ? Are they located in the very localized 
3d orbital of the transition metal, or in the extended 4s 
or even on the oxygen ligands ? In this paper we claim 
that resonant X-ray diffraction is able to answer quanti- 
tatively these questions. In our study, we have found a 
small but effective ordering occurring mainly in the 3d 
transition metal orbitals. We think that the same trend 
occurs in many compounds exhibiting CO which range 
in diversity and physical properties from colossal mag- 
netoresistance manganites, to the astonishingly complex 
Verwey Q transition in magnetite, to the CO inducing 
Spin-Peierls phase transition in NaV^Os ||. 

X-ray absorption spectroscopies (XAS) and resonant 
X-ray scattering are extensively used techniques to ex- 
plore and elucidate the electronic properties of materi- 
als 0, ■ I n these spectroscopies the incoming photons 
are absorbed by promoting a core electron to an empty 
intermediate excited state /. When the electron sub- 
sequently decays to the same core-hole, it emits a sec- 
ond photon with the same energy than the incoming one. 
This process is virtual and preserves the coherence of the 
electromagnetic wave. Set in diffraction conditions the 
associated experimental technique is known as diffrac- 
tion anomalous fine structure spectroscopy or resonant 



diffraction. When the absorption process is considered 
alone, and therefore loosing the coherence properties of 
the X-ray field, the experimental technique is called XAS 
or X ray absorption near edge spectroscopy (XANES) 
when considering the structure close to the edge. In this 
low energy range a higher sensitivity to the electronic 
configuration around the photoabsorber is expected be- 
cause the photoelectron probes allowed electronic states 
very close to the Fermi level as required by the selection 
rules resulting from the Fermi golden rule. Indeed by 
playing with the polarization of the incoming (and even- 
tually outgoing) photon beam, with the different chan- 
nels of transition (dipole, quadrupole) and with the ori- 
entation of the crystal with respect to the photon polar- 
ization direction, the probability of transition becomes 
highly angle dependent. This phenomena is illustrated 
by a tensor form of the scattering structure factors, and 
the analysis of each term of the tensors gives access to 
specific projections of the electronic states. Although 
complete calculations are performed we will restrict to 
the dipole approximation and the atomic scattering fac- 
tor (ASF) reads, 
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a and [3 are the index of the coordinates in the orthogonal 
basis. e l and e° stand for the incoming and outgoing 
polarizations and D a p is given by: 
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where Ef,E g and hu> are respectively the energies of the 
intermediate state, the core state and the photon; m e is 
the electron mass and T(Ef) includes the contribution 
of the inverse lifetime of all the states involved in the 
resonant process. The summation over / is limited to 
the unoccupied states. In diffraction, the intensity of the 
reflections is proportional to the modulus squared of a 
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weighted sum of the ASF: 
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where R a stands for the atom a position and Q is the 
diffraction vector. The atomic scattering amplitude is 
the sum of the Thomson factor fo a and the ASF. 

Of particular importance are the forbidden Bragg 
peaks that are extinct off resonance conditions because 
of the presence of glide planes or screw axis, if allowed 
by the space group of the given material. For these 
reflections, the diffracted intensity is related to the 
difference between ASF of the resonant (or anomalous) 
atoms. Therefore, these reflections are supposed to 
be much more sensitive to slight modifications of the 
crystallographic as well as the local electronic structure. 
In the absence of new superlattice peaks originating 
from a phase transition (typically ferrodistortive phase 
transitions), these forbidden Bragg peaks might yield 
relevant informations on the corresponding CO. 

This phenomenology has been extensively used in the 
study of the interplay between orbital, spin and charge 
degrees of freedom in transition metal oxides 
particular it is specially suited to study charge 
orbital orderings 0, 0, 0|. Nevertheless apparently 
contradictory results have led to opposite conclusions 
on the presence or not of real CO. Recent resonant 
diffraction experiments at forbidden reflections in Fe3 04 
have yielded no change in the resonant signal below 
and above the ordering temperature and it has been 
concluded that no CO occurs at the Verwey transition 
|l4| . Conversely, resonant scatte ring experiments in 
NaV205 0| and La 1 / 2 Sr 1 / 2 Mn03 [la] have shown that 
the energy line shape of some reflections retains close 
similarities with the derivative of the resonant X-ray 
scattering factor, a fact which has been interpreted as 
a signature for the presence of two different V (or Mn) 
charge states ordered in the lattice. In this letter we shall 
show that all these experimental evidences can actually 
be interpreted within a unique framework. In order to 
accomplish this task, we shall proceed in a quantitative 
analysis of the resonant diffraction experiments with the 
help of ab initio simulations. It will then be possible to 
evaluate the influence of the different parameters and 
estimate the effective charge redistribution resulting 
from the CO transition. 

To illustrate our analysis we shall discuss very recent 
resonant X-ray diffraction results in NaA^Os and all the 
experimental details can be found in ref. |l5j. Among 
the numerous diffraction peaks measured as a function 
of the incident photon energy and under different photon 
polarization conditions, we select three of them to illus- 
trate our results. The first one is the (1, 0, 0) c forbidden 



reflection which does not exhibit any change through the 
transition. The incoming photon polarization is along 
c, the outgoing one is in the scattering plane (a — tt 
polarization condition). The two other reflections under 
consideration are superlattice peaks, (15/2, 1/2, l/4)f, jC , 
for which the incoming and outgoing photon polariza- 
tions were kept parallel (a — a condition). According 
to the index, the polarization directions are along either 
the b— or c— axis. The superstructure peaks appear due 
to the phase transition and the energy dependence of 
the scattered intensity exhibits peculiar features that we 
examine in detail below. Ab initio monoelectronic calcu- 
lations are performed using the FDMNES code. 0, 0] 
Starting from the atomic positions and an associated 
electronic density, this code solves the Schrodinger or 
Dyson equations to determine the intermediate states 
probed by the photoelectron. The resulting structure 
factors are evaluated by taking into account the different 
cxcitation-dccxcitation channels (dipole-dipole, dipolc- 
quadrupole and quadrupole-quadrupole). Diffracted 
intensities are calculated with the crystallographic 
structures corresponding to the high (Pmmn |l9|) and 
low (All 2 ,20]) temperature phases with and without 
CO. Because the purpose of the present paper is to 
fit the electronic structure, it is important to define 
easy to handle parameters. We choose the valence 
orbital occupancy rates (VOOR). The radial shape 
of the atomic orbitals remains fixed. For any VOOR 
set, the total electronic density is obtained by simple 
superposition of the electronic orbital densities weighted 
by the VOOR. First-guess VOOR values corresponding 
to the high temperature (HT) phase are obtained by a 
fit on the electron density calculated by a LAPW \2l\ 
band structure procedure. These orbital occupancies, 
which reflect the location of the charge in space, are 
further refined in order to accurately reproduce the 
experimental curves in the low temperature phase. 
The vanadium and oxygen VOOR are then the only 
adjustable parameters of our calculations. The HT 
starting values are respectively 2.3 and 1.7 electrons 
for the 3d and 4s V atoms and 4.6 electrons for the 
2p O atoms. If one uses the ionic radii as a definition 
of the atomic spheres, the 4s charge (average radius 
is 1.5 A, to be compared with the 3d average radius, 
0.6A) lies inside the oxygen atoms. Thus, the V formal 
charge is 2.7+ in this context. The sodium atom has no 
electron on its 3s valence orbital and remains 1+ in all 
the study. With this VOOR set, a good agreement is 
obtained at HT between the calculated (100) reflection, 
the polarized XANES line shape and the corresponding 
experimental spectra. To get such agreement, the finite 
difference procedure is used to solve the Schrodinger 
equation. This technique permits to avoid the muffin-tin 
approximation on the potential shape. That point is 
absolutely necessary in such distorted systems. 
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We now discuss the low temperature (LT) phase and 
the starting point is the crystallographic structure as 
determined in ref. [l9|. The most important result is 
that an explicit CO is absolutely necessary to get a 
satisfactory agreement with the experiment (Fig. 
The fit of the experimental data yields two equivalent 
solutions : (a) The charge exchange is located at the 
4s vanadium and amounts 5 = ± 0.2 electron, and 
(b) it is located at the 3d vanadium and amounts only 
8 = ± 0.04 electron. Thus the CO takes place either 
in the intermediate region between the oxygen and the 
vanadium where an atomic assignment of the electron is 
merely formal (see Fig. EJ or closer to the atom centers 
with a smaller magnitude. In any case, the bridging 
oxygen located in the rung of the ladder, being the 
atom that experiences the largest displacement at the 
phase transition (0.05A towards the V 5+ , that is 5 
to 10 times more than the other atoms), is suspected 
to participate to the electronic rearrangement as well. 
Thus, the electronic reorganization we propose does not 
concern isolated atoms but all the rung of the ladder, 
with the center of gravity of the electronic density 
displaced in the opposite direction than that of the 
bridging oxygen, i.e., toward the V 4+ atom. Taking 
into account the experimental resolution, the choice 
among the two proposed solutions cannot be safely 
done. Nevertheless the second one with a CO on the 3d 
orbital is in agreement with that obtained by Suaud and 
Lepetit on the basis of ab initio calculations using 
large configurations interaction method on embedded 
fragments. These authors found a ± 0.04 CO on the 
3c?, in excellent agreement with the present study. 
They have also shown that the modification of the 
magnetic properties at the phase transition is due to the 
redistribution between the different d— orbitals, keeping 
the total d charge only slightly modified. We want to 
end this paragraph by underlying the extreme sensitivity 
of resonant diffraction to observe such a small amount 
of charge disproportionation. 

A second important result, with straightforward im- 
plications in other CO compounds, is that the CO has 
very little influence on the (1, 0, 0) forbidden reflection 
whereas it is at the origin of the derivative line shape 
for the non integer reflections. A charge disporportiona- 
tion in the vicinity of the absorbing atom greatly affects 
localized orbitals, and thus both the Is and the 3d en- 
ergy levels will undergo an analogous shift. However the 
4p states, less localized, are hardly affected j^. There- 
fore the dipole-dipole contribution to the K— edge of an 
atom experiencing a CO will be shifted towards higher or 
lower energies as compared to that found in the HT phase 
whereas the quadrupole pre-edges are little shifted dur- 
ing the oxidation processes. In a more general way, one 
has to proceed very cautiously on applying the derivative 
effect concept to L— and M— edges if the CO is expected 
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FIG. 1: Experimental and calculated resonant diffracted peak 
spectra in NaV^Os. From bottom to top the (1, 0, 0) a — n, 
(15/2, 1/2, l/4) c and (15/2, 1/2, l/4) b peaks (both measured 
in the a — a configuration) and the XANES cross section are 
shown. For each reflection the spectra corresponding to the 
experiment (dotted line) and the calculation in the low tem- 
perature phase with charge ordering (continuous line) and 
without charge ordering (dash line) are displayed. For the 
(1, 0, 0) reflection the high temperature structure calculations 
are not shown as they are practically identical to the LT cal- 
culation without charge ordering. The charge order model 
yields very few modifications of the signal through the tran- 
sition for the first peak. The derivative effect for the two last 
ones can not be accounted for by a lattice distortion alone, 
and is best reproduced when different charges are introduced 
at convenient V-sites. The XANES spectra for the polariza- 
tion e along the c— and axis are also given. The different 
calculations at HT and LT with or without charge ordering 
yields practically the same results. 



to take place in d— and /— final states, respectively. We 
have computed the effect of the associated geometrical 
distortion onto the ASF and we have concluded that it 
remains negligible (± 0.1 eV) as compared to the ob- 
served energy shift and to the effect of the charge dispro- 
portionation (SE = ±0.8 eV). A first order Taylor ex- 
pansion of the scattering factor is justified by the small 
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FIG. 2: A sketch of the NaV20s structure focussed around 
the V 4+ -0-V 5+ rungs. All oxygen atoms lie in the same plane, 
whereas vanadium rungs lie alternatively below and above the 
oxygen plane. The vanadium (full circles) are surrounded by 
five oxygen atoms, four shown in white circles on the figure 
and a fifth one in the apical direction (not shown) . The rings 
show the region were the electronic density is affected by the 
CO, the difference in density around the vanadium atoms be- 
ing illustrated by the different grey tones. 



value of SE. It reads as f+ 5 « / + df/dE ■ SE and 
f- 5 w / - df/dE ■ SE where the /, f +s and /"* are the 
ASF tensors for the vanadium atoms with the HT charge, 
and the +<5 and —5 charge disproportionation, respec- 
tively. The derivative effect condition is best achieved 
for diffraction wave vectors with direction and amplitude 
such that the putative disproportionated vanadium ASF 
subtract. Conversely, for the forbidden reflections -as the 
(1,0, 0)- the derivative term remains negligible in front 
of the combination of the ASF tensors weighted by its 
phase factors, and the sum of phase factors of V-atoms 
of a given charge state almost cancels out. 

By using the multiple scattering theory (MST) we have 
checked that, except for the energy shift, the ASF are 
very little modified through the transition. Therefore 
our assumption above is plainly justified and we use en- 
ergy shifted HT ASF in the LT structure. It can also be 
noticed that the XANES spectra show little sensitivity 
to the charge occupancy. The reason is exactly the same 
as the non-sensitivity of the (1, 0, 0) reflection: XANES 
adds the contributions of the different atoms, the small 
shifts due to the CO remain smeared within the experi- 
mental resolution. 

We have shown that a charge-ordering transition model 
quantitatively explains the resonant scattering experi- 
ments in NaV205. The method presented here has the 
advantage to allow for a direct evaluation of the actual 
ordered charge disproportionation in the 3d-transition 
metal oxides. The charge ordering concerns a small frac- 
tion of an electron mainly the d— orbitals even if a con- 



tribution of the s— orbitals cannot be excluded from the 
present study. This conclusion, although heavily sus- 
pected from the results of other works, points to recon- 
sider electronic interactions in strongly correlated elec- 
tron systems and to the role of small lattice distortions 
in accommodating extra charges. Finally with respect 
to compounds such as manganites and magnetite, where 
resonant diffraction experiments have yielded no sign of 
CO, we suspect that the studied diffraction peaks are 
hardly sensitive to an actual CO. Our results in NaV205 
give evidence that spectroscopic effects associated to a 
very small CO are not at all negligible and that they 
should be observable in resonant diffraction experiments 
if convenient meaningful peaks are measured and the re- 
sults contrasted with the help of appropriate simulations 
codes. 
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